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The analytical method of solving the boundary problems for a system 
O ■ of equations describing the behaviour of electrons and an electric field in 

C ■ the Maxwell plasma half-space is developed. Here the diffusion reflection of 

electrons from the surface as a boundary condition is supposed. The exact 
(N ; solution for boundary value problem of the Vlasov - Maxwell equations 

through the eigenfunctions of the discrete spectrum and the continuous spect- 
rum is obtained. The method of decomposition by eigenfunctions of discrete 
and continuous spectrum is used. The exact expression for the surface impe- 
dance also is obtained. 
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(N : PACS numbers:52.35 - g; 52.20 - g; 52.25 - b 

> 
oo 

g ■ 1. Introduction. Statement of problem. The skin effect is the plasma's 

response to an external constant-amplitude variable electromagnetic fleld 

o ■ 

0\ ; that is tangential to the surface [Ij. This problem represents essential interest 

O ■ 

with the focal point being on the computation of surface impedance. 
^ ■ In 1^ showed that electrons can transport the plasma current away from 

5^ i the skin layer due their thermal motion over distances of order vt/co. In 

[S] this numerical study provides through speciflc numerical comparisons a 
useful perspective on the limitations and capabilities of the various Krook 
collisional models. This information should be useful to future studies that 
try to flnd a compromise in handling the yet unresolved problem of hoe to 
unify the kinetic response and collisions in wave phenomena. 

The paper further develops the analytical method of solving boundary 
problem for a system of equations describing behaviour of electrons and 
electric fleld in Maxwell plasma. The method is based on the idea of the 
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decomposition of the solution the problem by singular eigenfunctions of the 
corresponding characteristic equations. 

We find the analytical solution to the boundary problem of the skin effect 
theory for electron plasma that fills the half-space. The analytical solution 
has a form of the sum of an integral and two (or one) exponentially decreasing 
particular solutions of the initial system. Depending on index the problem 
one of the particular solutions disappears. 

We consider a collision plasma with behaviour is described by modelling 
kinetic equation with integral of collisions in the form BGK-model (BGK= 
Bhatnagar, Gross, Krook). 

Let's Maxwell plasma fills the half-space x > 0. Here x is the orthogonal 
coordinate to the plasma boundary. Let's the external electric field has only 
y component. Then the self-consistent electric field inside in plasma also has 
only y component Ey{x, t) = E{x)e~^'^^. We consider the kinetic equation for 
the electron distribution function: 

| + .,| + eS(.)e-"g = .(/o-/(,..v)). (1.1) 

In equation (LI) i/ is the frequency of electron collisions with ions, e is the 
charge of electron, /o(v) is the equilibrium Maxwell distribution function, p 
is the momentum of electron, 

/o(v) = nf^V^'exp(-/3V), = 



TV J ^ ^ 2kBT' 

Here m is the mass of electron, is the Boltzmann constant, T is the 
temperature of plasma, v is the velocity of the electron, n is the concentration 
of electrons, c is the speed of light. 

The electric field E{x) satisfies to the equation: 

E"{x) + -^E{x) = / Vyfit, X, v) d\. (L2) 

We assume that intensity of an electric field is such that linear approxima- 
tion is valid. Then distribution function can be presented in the form: 

/ = /o (1 + C'y exp(-ia;t)/i(x, /i)) , 
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where C = is the dimensionless velocity of electron, ji = Cx- Let 

/ = vtt is the mean free path of electrons, vt = vt is the thermal 

electron velocity, r = 1/z/. We introduce the dimensionless parameters and 
the electric field: 

X . . y/2e _,, , 
t-^ = jyt, xi = -, e{xi) = — -==E{xi). 

Later we substitute xi for x. The substitution produces the following form 
of the kinetic equation (1.1) and the equation on a field with the displacement 
current (1.2): 

dh 

li— + zoh{x,ii) = e{x), zo = l-icjT, (1.3) 

00 

e"ix) + Qh{x) = -i^ I exp(-//'') h{x, //') d^jl , Q = —, (1.4) 

J c 

— DO 

where 6 = , 6 is the classical depth of the skin layer, ctq = , ctq is 

zTTuao mu 

2/2 

the electric conductance, a = -z^, a is the anomaly parameter. 

We formulate the boundary conditions for the distribution function of the 
electron in case of the diffusion electron reflection from the surface: 

/i(0,//) = 0, 0<//<+oo. (1.5) 

We use the condition that function h{x,iJ.) vanishes far from the surface: 

/i(+oo, yLi) = 0, — oo < /i < +00, (1.6) 

and conditions for electric field on the interface and far from it: 

e'(0) = e/, e(+oo) = 0, (1.7) 

where is the given value of the gradient of the electric field on the plasma 
interface. 

So, the skin effect problem is formulated completely. We seek solution of 
system of the equations (1.3) and (1.4) in this problem, that satisfy to the 
boundary conditions (1.5) - (1.7). 
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2. Eigenfunctions and eigenvalues. The separation of variables |[2]-[1] 

X 

h^{x, /i) = exp(-zo-)^(T/, /i), (2.1) 

X 

er^{x) = exp{-zo-)E{T]), (2.2) 

where r] is the spectral parameter (generally, complex- valued) , reduces system 
of the equation (1.3) and (1.4) to the characteristic system: 

(r/-^)$(T/,/i) = ^E(T/), (2.3) 
[4 + QY]E{r,) = -'^n(v), (2.4) 



where 

n{r]) = / e~^^^{r], fi)dfi. 



-oo 



From equations (2.3) and (2.4) for 77 G (—00, 00), we find the eigenfunctions 
of the continuous spectrum [2], [B]: 

$(r/, /i) = ^rj^e-^'P— + A(^)(5(r/ - /i), (2.6) 



E{r^) = ^r^'e-^\ a=-i-,. (2.7) 

In the equation (2.6) the symbol Px^^ denotes the distribution, i.e. the 
principal value of the integral of x~^, 6{x) is the Dirac delta function. 

Since the distribution and electric field functions decline as we go further 
from the boundary, we regard the continuous as positive real half- axis: < 
T] < +00. The eigen solution of the continuous spectrum hn{x,fi), e,y(a:) are 
decreasing from value since Re > 0. 

In the equations (2.6) and (2.7) there is the normalization condition: 
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J e ^'$(r/, ii)dii = 

—00 

X{z) is the dispersion function. 



1. 



X(,) = l+(Q)\^ + ^f!:l±. (2.8) 
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Consider the dispersion function in the following form: 



X{z) = l + bz — ap{z), 
here b = ^ = ^^, = / ^^EiziMi. 



6 = 



a = —I- 



The fact that A(/i) has a double pole at 2; = (X) follows from the asymptotic 
series in neighborhood of infinity: 

X / N N 9 a\ 3a 1 15a 1 .^.^x 

3. Mode of plasma in Maxwell plasma. Find out the structure 
of discrete spectrum, this spectrum consists of the zeros of the dispersion 
function. Each zero corresponds to the own decision of the discrete spectrum, 
named also mode of plasma. 

Take two line F^, parallel real axis and defending from it on distances, e > 
0. The value e shall choose so small that all zeroes to dispersion function lay 
outside of narrow band, concluded between direct F+ and F~. 

According to principle of the argument difi"erence between number of the 
zeroes and number pole to dispersion function is an incrementation of its 
logarithm: 

1 



N -P = 



2m 



d\n\{z). (3.1) 



'■r+ Ti 

In (3.1) each zero and pole are considered so much once as their multiplicity, 
the lines F+ and F~ are passed accordingly in positive and negative directions. 
Then the dispersion function in infinitely removed point has a pole of the 
second order, i.e. P = 2. 
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When £ — )• from (3.1) we get: 



N -2 = 



2m 



din 



a+(m) 



— oo 



Integral from (3.2) represented as: 



oo 



/ 



dhi 



A+(^ 

a-(m) 



A+(^) f X+{ii) 



(3.2) 



—00 — 00 

In the second integral we shall do change the variable: r — r, we have: 



A+(-r) = A-(r), A-(-r) = A+(r). 



Consequently, second integral is first. Really, 





oo 



A-(/^) 



— OO 



A-(-//) J 





oo oo 



A+(//) 



A-(//) 



Thereby, 



00 



M-2 = lfdin^ 

TTl J X-{ll) 





(3.3) 



On complex plane we consider line F : z = G{t), ^ t ^ +oo, where 
G{t) = ^-4^. Not difficult check that 



x-it) 



G(0) = 1, lim G{t) = 1. 



These equality mean that curve F is closed: it comes out of points z = 1 
and ends in this point. According to (3.3): 



N -2 = — 



ln|G(r)| +iargG(r) 



+0O 1 
TT 



arg G{r) 



+00 



Thence we get: 



N-2 = - 

TT 



- +00 

argG(0 =2x{G), 



(3.4) 



where x = >€{G) is the index of the problem, G{t) is the number of revolution 
the curve F comparatively begin coordinates, made in positive direction. 
From (3.4) we see 

Ar = 2 + i[argG(+oo)-argG'(0)l = 2 + i arg ^(+00), (3.5) 

TT L J TT 
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since arg G{0) = 0. 

argG(t) is the regular branch of the argument, fixed in zero condition: 
arg G{0) = 0, and determined in cut complex plane with cut along positive 
part of real axis. This branch complies with the main by importance of the 
argument. 

We present the dispersion function in form (the displacement current is 
absent): 

oo 

, ri f exp(— . , , 

X(z) = -zc - + ^ / r J r ^ -^cLj(z), 

Lc ^/^^ J ii - z J 

— 00 

where 

CXD 

cj{z) = 6 + ^ / , 6 = -. (3.6) 

Vtt 7 II- z c 

— oo 

According to (3.6) border values of the function lo(z) overhand and from 
below on real axis are: = cj(/i) i iy/i^ fi^ ex.Y>{~ 1^^) ■ 

Then, the G{ii) is equal: 



Where 



whence 



4c^ (cu + iiyf 
5 = 61 + 162 = 




TTVt 



^1 = 2 2~' ^2 = '^c = . 

Select in equality (3.7) for function G(//) real and imaginary part: Gi^ji) = 
Gi{ii) + iG2{lj)- We have 

[Ji-KAi)]' + ^i-9'(Ai) 



Gi(/^) = 



[Si - pif^W + [S2 - qif^W 
2q{,j) [5i-p{li)] 



On complex J-planes shall enter curve A: 

A: (5i=]5(^), 52 = ±q{ii), ^ /i ^ +00. 
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Entrails the curve A shall mark through A"*", and A~ is the external of 
curve A. The domains A"*" and A~ are unlimited(fig. 5.1). 
Possible prove that: 

1) if ^ G A+, that arg G(+oo) = 27r (the curve A does one turn in positive 
direction comparatively begin coordinates), 

2) if (5 G A~, that argG(+oo) = (the curve A does not cover begin 
coordinates) . 

0.6 

0.3 

0.0 
-0.3 
-0.6 

Fig.l. The domain A='= on (5-plane, S — {Si, 62). 
According to (3.5) we get 

1) if J G A+, then = 4 - dispersion function has four zeroes, and 

2) if (5 e A~, then N = 2 dispersion function has two zeroes. 

Go from (5-planes to planes (a;i,z^i), where uji = 'y/vc, = e/vc, and 
7 = uj/ujp, e = u/ujp. Find images the unlimited domain A^ and curve A 
under such image. 

Present the parametric equations crooked A in form 

Li : -cji + Sz/j = p(yw), -2>uJiUi + — = ±q{ii). (3.8) 

From the first equation from (3.8) find: 

i^i = ±^H-p{ii). (3.9) 
Substitute (3.9) in the second equation from (3.8). We Get: 

vl - SuiiSi^l - pill)) = ±q{ii)^?>ul-p{ii). 
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Fig.2. Unlimited domain Df on (cui, z/i)-plane. Border of the domain is the curve 

Li = dDf. 

Involve this equation in square: 

64iyt - 48iytpil^) + Sz/^ [Sp\,^) - q\,^)] + q\ii)piii) = 0. (3.10) 

Under each fixed ji G [0, +oo] the equation (3.10) has a single zero, which 
we mark through Y[fj,). Ensemble of all such zeroes forms is function z/i = 
y(//), O^fj,^ +00. 

Hence we find the first parametric equation of curve Li: 

e = VcY{lJ,), O^yLi^+OO. (3.11) 

Substituting (3.11) in the second from equations (3.8), find: 

oji = v^3y2Qiy^r^(/I), < // < +00. (3.12) 

According to (3.11) and (3.12) we get parametric equations curve Li on 
planes {uJi, ui): 

Li : uji = ^/3Y'^{|J,) -p[ijl), i^i = Y[ii), ^ /j, ^ +oo. 

Go from plane {cui, vi) to planes parameter (7,^). Present the parametric 
equations crooked A in form: 

^3 

-7^ + 3£^ = vIp{ii), -37£ H = ±Vcq{ii). (3.13) 

7 
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where 

7 



Up LUp 



From the first equation from (3.13) find: 

-f = ±^Se^-v^pifj,). (3.14) 
Substitute (3.14) in the second equation from (3.13). We Get: 



Involve this equation in square: 

6 /P\4 / P \ 2 



64(^) -48(^) K)^) + 3(^) [?.p\^i)-q\^i)\+q\^i)p{^i) = 0. (3.15) 



Under each fixed G [0, +oo] the equation (3.15) has a single zero, which 
we mark through Y^ji). Ensemble of all such zeroes forms is function 

Consequently, we have found one parametric equation of curve L, being 
image of curve A: 

e = VcY{^i), 0^)Li^+(X). (3.16) 
Substituting (3.16) in the second from equations (3.13), we find: 



-i = Vc^/2>Y\^l)-p{^l), 0^//^+oo. (3.17) 

According (3.16) and (3.17) on (7, £)-plane we have the parametric equati- 
ons of curve L{vc)- 

L{vc) : 7 = Vc\/'iY'^{ii) - pill), e = VcY{ii), ^ II ^ ^00. 

We find the general solution to system (1.3), (1.4) in the form of an 
expansion in the eigenf unctions of the discrete and continuous spectrum; 
this solution automatically satisfies the boundary conditions at infinity: 

k=0 



h{x, )w) = ^ V exp - ) + 



00 



+ / exp(^-^)AinMn,i,)dri, (2.10) 



7] 




11 




0.0006 0.0012 0.0018 



Fig. 3. The unlimited domain D"^ on the (7, £:)-plane. For curves 1,2,3 
= 0.0006 (T = 1000°/^), 0.001 (T = 3000°/s:), 0.013 (T = 5000°/^) respectively. 



fc=0 

00 



+^ ei,Y>{-rf-^-^)rfA{r))dr]. (2.11) 


Here Ak{k = 0, 1) are the unknown coefficients corresponding to the 
discrete spectrum, = 0, if ^ G and A{r]) is an unknown function 
called the coefficient of the continuous spectrum Re (z/rjk) > (A; = 0, 1), 
Re zo = 1. 

4. The index of the problem equals zero. We consider the case when 
the problem index is equal to zero. Substituting expressions (2.10) and (2.11) 
into (1.5), we obtain the singular integral equation with the Cauchy kernel: 



00 



a f^^Mr])A{ij)_^^^^^^^^^^^^^^^^^^^^ 0<//<oo, (3.1) 



7T J ^ ~ 





where 
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We define the auxifiary function in the complex plane 



1 f exp(— 77^)74(77) 



TT 



dij. 



Then the equation (3.1) is possible to transform to the Riemann 
boundary value problem: 



Hilbert 



X+ifi) N+{fi) + cp{ii) = X-{ii) N^{ii) + ip{ii) , 0</i<oo. (3.2) 

We consider the Riemann — Hilbert boundary value problem on the half- 
axis: 



X-ifi) 



= G{ii), < /2 < +00. 



where 



G(/i) = 



A+(^) 



A-(/i)' 

G{ii) is the coefficient of the boundary value problem. 

The index in this problem is >£{G) = 0, then the solution of problem (3.2) 
has the form [4j: 



X{z) = expV{z), V{z) = ^^ J 



1 flnG{r)dr 



T — Z 



(3.3) 



We use (3.3) to reduce problem (3.2) for determining an analytic function 
N{z) from its zero jump problem on the cut: 

X+{ii) \n+{ii) + ip{ii)] = X-{ii) \n-{ii) + Lfiii)] , < /i < 00. (3.4) 



A general solution of (3.4) is given by the formula 

Co 



[z - m)xiz) 



v/tF z-r]o 

where Gq is an arbitrary constant. 

Let's eliminate poles at the decision (3.5), we have: 

^0 



(3.5) 



Gn = 



exp(-r/o)r/JX(r/o 



(3.6) 
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Employing the Sokhotsky formula for difference from the boundary values 
of function N{z), we find the coefficient of the continuous spectrum: 



^^irf exp(— 77^) ^4(77) = 



1 1 



77 - ryo X-{r]) 
From the definition of the auxiliary function ^{z) we have: 



00 

A^(0) = ^ / 77^exp(-772)A(77)(i77. (3.7) 


Considering equality (3.7), we write (2.11) in form: 

-^^ + ^(0) = -. 

r]oX{r]o) azo 
From general solution (3.5) the constant Cq is found: 

Co = -^X(O). (3.8) 

CLZq 

From comparison (3.6) and (3.8) we get the unknown coefficients corres- 
ponding to the discrete spectrum: 

° azoX{r]o)r]^ e^p{-r]^)' 
The derivative of expression (2.11) leads us to equality 



e'{x) Co zqx 
exp( ) 



00 



/ exp(— 77^ )r)A(r])dr]. 

/tt J 77 



az^ rilX[T]Q) 770 ^/7^ 



From here we find the electric field derivative on plasma boundary x = 0: 

ex.p{—r]^)r]A[r])dr]. (3.9) 



00 



e^(0) _ Co _ _i_ I 
az^ r}oX{if]o) 





For calculation of integral from the expression (3.9) we shall use a derivative 
of auxiliary function 

00 
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Using the expression (3.10) we find the derivative in zero: 

1 



CX) 



N'{0) = 



T] exp(— ?7 )A{r]) drj. 



(3.11) 



Using expression for coefficient of the continuous spectrum by (3.11) we 



find: 



whence 



Cn 



2-r/o)HX(z) X(t/o 



^'(0) = -§( 1 



7/2VX(0) X(r/o)^ 7/0^2(0) 



Co X'{z) 
Co X'{0) 



(3.12) 



Substituting (3.12) into (3.9) we obtain the equation: 



e'(0) 



azpCp 



1 



X(0) 



= zo 



X(0) 



1 

^0. 



(3.13) 



In [2] the following formula for the impedance is given: 

_ 4:mul e(0) 



(3.14) 



c^zo e'(0) 

Substituting (3.13) into (3.14), we obtain the exact expression for the 
impedance: 

'x'{o) _ r 

^0. 



z = 



4:7TiUjl 



c2 Zo 



X{0) 



-1 



(3.15) 



The equation (3.15) is the exact expression for calculation of the surface 
impedance. This equation expresses value of the impedance in terms of the 
function X{z) and zeros of the dispersion function of the problem. 

5. The index of the problem equals one. Now we consider case when 
the problem index is equal one, i.e. x(G) = 1. 

Substituting expressions (2.10) and (2.11) into (1.5), we obtain the singular 
integral equation with the Cauchy kernel: 



oo 



a(p{fi) + 



where 



a f T] exp{—r] )A{r]) 



T]- 



df] + X{ii)A{fi) = 0, < ^ < oo, (4.1) 



fc=0 



r]k- 
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We define the auxifiary function in the complex plane: 

00 

1 -tn'^M-nMin) 

^J-K J 7]- Z 



for which boundary values from above and from below on the valid axis 
formulas of Sokhotskii are carried out: 

Ar+(//) - Ar-(//) = 2V7ri//^exp(-//2)A(//). 

Using the boundary values N^ifi) and A^(//), we pass from the singular 
equation (4.1) to the Riemann — Hilbert boundary value problem: 

X+ifi)[N+{f,) + ^(/i)] = X-ifi)[N-if,) + ^(/i)], < // < oo. (4.2) 

Let's solve a corresponding Riemann — Hilbert boundary value problem: 
X+(/i) 



X-(/i) 

where 



= G(//), < // < +00, (4.3) 



Since problem index is equal one x(G) = 1 as the solution of (4.3) we 
take function 

oo 

X{z)=^-e^pV{z), Viz) = ^ / '"^<^)-^" rfr. (4.4) 



Using (4.4) we transform the boundary value problem (4.2) to determining 
an analytic function from its zero jump problem on the cut: 

X+(//)[7V+(//) + ^(//)] = X-(//)[7V-(//) + ^(//)], < // < oo. (4.5) 

The general solution of (4.5) is following 

where Co, Ci is the arbitrary values. 

Since A^(oo) = 0, from here we receive that Co + Ci = 0. 
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Let's eliminate poles at the decision (4.6), we get 



Cn = 



4eM-r]l)Xivo), Ci = -^T^leM-lDXM- (4-7) 



TT 



Substituting the general solution (4.6) in the Sokhotskii formula for (4.6) 
we find unknown function called the coefficient of the continuous spectrum: 

1 



77^ exp(— 77^)^4(77) = 



1 



1 



1 



Using (4.7) we reduce (2.10) to the form: 

00 

Co Ci 



(4.8) 



If 1 
— + ^ I rf ex^{-if)A{ri)dT] = — , (4.9) 
1) VTTj azo 



From equation (4.9) we will find last integral using definition of function 
Niz): 



N(0) = ^ I r]'^exp(-r]'^)A(r])dr]. 



Using the common decision (4.6), we found A^(0): 
iV(0) = Co ^ ^ ^ 



+ 



.77oX(77o) 771^(^1) mXiO) 77iX(0) 
We substitute this expression into the equation (4.9): 

X(0) 



azo 



(--- 

V7?i 77n 



(4.10) 



Differentiating expansion for electric field (2.11) and using (4.7), where 
a; = we found: 

CXD 

e'(0) = - - ^ [r^exp{-r^')Ain)dfj. (4.11) 





Integral from (4.11) is the auxiliary function derivative A^'(O): 

00 

iV'(O) = A=[ r]e-xj^{-if)A{r])dr]. 
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The common solution (4.6) we present in form: 



Niz) = 



Cn 



z -rio iX{z) a/tF 



+ 



1 



z-iji 

Derivative of this expression is equal: 



N\z) = 



1 



{z - IJoY 

Co 



1 



Co X'jz) 
X{z) X(77o)J z-77oX2(z) 



1 



1 



Hence, we obtain: 



+ 



Co X\z) 
z — rji X'^{z) 



N'iO) = -\ 
To 



1 



1 



+ 



Co 

7)1 



1 



X(0) X(/7o) 
1 1 



+ 



Co X'{{)) 

/?0X2(0) 

Co X'(0) 



X(0) X{r]i)\ 77iX2(0)- 
Using this equality for the electric field derivative we have: 



e'(0) = 



X(0) 



Therefore according to (4.10) we receive that 



e'(0) = zo 



X'(0) 



1 



1 



Hence, expression for the impedance is equal: 



Z = 



c^zo 



X'{0) 1 1 

'x{oj~^o~m 



-1 



This equation expresses value of the impedance in terms of factor function 
X{z) and zeros of dispersive function of the problem. 



CONCLUSIONS 
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A closed form solution of a system of two equations Boltzmann — Vlasov 
and Maxwell arising in the skin effect problem for a Maxwillian plasmas is 
presented. The kinetic Boltzmann — Vlasov equation with a r-model collision 
operator is considered. Case's [7] method and Riemann — Hilbert boundary 

value problem [8] with the coefficient G{fi) = -^-Jr-l, where X{z) is the 
dispersion function of the problem. The discrete fashion is obtained. Consider 
the domain that if the frequency lies in this domain, there are four 
discrete solution, and if the frequency is outside this domain, there are two 
discrete solutions. The exact solution of the initial boundary value problem 
with diffusion scattering of electrons from plasma boundary is constructed. 
The exact formula for calculation of the surface impedance is obtained. 
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